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Altitude 


Sao  Luis,  Jicamarca  and  Christmas  Island  plumes 
caracterization  in  function  of  solar  flux 

The  following  plume  parameters  were  analyzed  for  Sao  Luis,  Jicamarca  and  Christmans  Island: 


Range- Tims- Intensity  map,  01  December  2002 
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Hi  -  Bottom  type  onset  altitude 
Hp  -  Plume  onset  altitude 
Hpk-  Peak  plume  altitude 
Ti-UT  Time  for  Hi 
Tp-UT  Time  for  HP 
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Plumes  caracterization  in  function  of  solar  flux 

SAO  LUIS  2001-2003  (EQUINOX) 
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Onset  altitude  of  bottom-type  ,  plume  and  plume  peak  increase  almost  linearly 
with  increasing  solar  flux.  Onset  time  decrease  with  solar  flux. 


Plumes  caracterization  in  function  of  solar  flux 


SAO  LUIS  2001-2008  (NOV-FEB  Summer  solstice) 
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Onset  altitude  of  bottom-type  ,  plume  and  plume  peak  increase  almost  linearly 
with  increasing  solar  flux.  Onset  time  presents  almost  no  variation  with  solar  flux. 


Plumes  caracterization  in  function  of  solar  flux 


Bottom-type,  plume  onset  and  peak  altitude  variation  with  solar  flux  for  the  3  stations 


EQUINOX  |  NOVEMBER- FEBRUARY 


Bottom-type  and  plume  onset  altitude  present  similar  behavior  for  the  3  stations,  while  the 
plume  altitude  peak  presents  smaller  inclination  for  Christmas  Island. 
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Plumes  caracterization  in  function  of  solar  flux 


Bottom-Type  and  plume  onset  time  with  solar  flux  for  the  3  stations 
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Bottom-type  and  plume  onset  time  presents  similar  behavior  for  Sao  Luis  and  Jicamarca 
however  larger  variation  (onset  time  earlier  for  higher  solar  fluxes)  for  Christmas  Island. 


Plumes  caracterization 


Ionospheric  irregularity  zonal  velocity  calculation  using  VHF  radar  interferometry  (in 

collaboration  with  Fabiano  Rodrigues  from  UTD) 


Sao  Luis  VHF  radar  antenna  sets  experimental  setup 
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Plumes  caracterization 
Sao  Luis  VHF  radar  experimental  setup 


150-km  echoes*  F-region* 


TX  power 

4  kW 

2x4  kW 

Code  length 

28  bauds 

28  bauds 

Baud  length 

1  km 

2.5  km 

IPP 

600  km 

1,400  km 

Initial  sampling  height 

90  km 

200  km 

Max  sampling  height 

210  km 

875  km 

Doppler  velocities 

±  40  m/s 

±  268  m/s 

♦Typical  configuration 
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Plumes  caracterization 
Results 


03- Dec-2010 


07-Dec-2010 


Sao  LUIS  -  03-D6C-2010  Sao  UllS  -  07-Dec-20i0 


Local  Time  Local  Time 
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Plumes  caracterization 
Results 


Relatively  small  height  variation  03-Dec-2010 
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Plumes  caracterization 
Results 


Slower  drifts  at  lowest  heights.  07-Dec-2010 


Conclusions  for  zonal  drift  calculations 


•  Results  confirm  the  possibility  of  estimating  relative  zonal  velocities. 

—  Regions  of  westward  drifts  (bottomside  F-region). 

—  Determination  of  the  height  of  the  shear  node  (as  high  as  450  km). 

•  The  turbulent  nature  of  ESF,  however,  makes  it  less  than  straightforward 
to  accurately  determine  the  absolute  mean  bulk  flow. 

•  Future  work  includes: 

—  Long-term  (several  minutes)  tracking  of  irregularity  dusters  to  derive  mean  zonal 
drifts. 

—  More  comprehensive  comparison,  if  possible,  with  plasma  bubble  velocities  and 
winds. 


GPS  signal  analysis 


S4  occurrences  (in  collaboration  with  Dr.  Alison  Moraes  from  IAE/CTA) 


Distribution  of  the  S4  indices  of  the  observations  available  between  Dec  14th  2001  and  Jan  14th 
2002,  as  a  function  of  local  time. 


GPS  signal  analysis 

Decorrelation  time 

The  autocorrelation  function  of  the  normalized  signal  amplitude  scintillation  is  given  by 

<  ,  ,  E[(R(t)-z)(R(t  +  r)-z)] 

Ar(t)  = - 5 - 

aR 

The  x0  value  is  defined  as  the  time  lag  at  which  the  autocorrelation  function  falls  off  by  e' 
from  its  maximum  (zero  lag)  value: 


While  the  S4  index  is  an  indicator  of  the  depth  (or  magnitude)  of  amplitude  fadings,  the 
decorrelation  time  (x0)  is  an  indicator  of  rapidity  of  the  fadings. 


C/NQ  (dB-Hz)  C/NQ  (dB-Hz) 


GPS  signal  analysis 

The  variability  of  amplitude  scintillation  patterns 

xq=0.94  (seconds)  xQ=0.68  (seconds) 


Examples  illustrating  the  variability  in  the  decorrelation  time  made  during  the 
campaign  of  observations  used  in  this  study.  In  all  cases  the  time  series  of  measured 
signal  amplitude  have  approximately  the  same  S4  (~0.9),  but  very  distinct  x0  values 


GPS  signal  analysis 

Example  of  observations 


Local  Time  l_oca|  Tjme 


Sao  Jose  dos  Campos  on  December  14th,  2001.  Panel  (a)  C/NO  of  the  LI  signal  received 
from  satellite  PRN  28.  (b)  S4  index,  (c)  azimuth  and  elevation  angle  of  the  PRN  28.  (d) 
decorrelation  time  (x0)  (is  not  computed  for  S4<0.1).  This  example  illustrates  that  the 
scintillation  intensity  decreases  and  the  decorrelation  time  increases  as  time  progresses. 


E[tq]  (seconds) 


1.1 


GPS  signal  analysis 

S4  vs  x0 


•  Linear  relationship  between  S4  and  xO; 

•  Variability/spread  of  xO  values  as  well  as  the  mean  value  of  xO  tend  to 
decrease  as  scintillation  intensity  (S4)  increases. 

•  Decorrelation  times  decreases  as  scintillation  becomes  stronger,  and 


18 


vice  versa. 


GPS  signal  analysis 

Amplitude  distributions 

S4=0.6  S4=0.7 


—  Scintillation  data 

—  Nakagami-m 

—  Rice 


ot-p  distribution  (Yacoub,  2007).  This  model  assumes  the  signal  as  a  composition  of 
many  clusters  of  multipath  waves  instead  of  just  one.  The  result  of  such  an  assumption 
is  a  more  comprehensive  characterization  of  electromagnetic  scattering  phenomena. 
Our  results  show  that  the  Nakagami-m  PDF  performed  better  than  the  Rice  PDF 
a-p  distribution,  however,  outperformed  the  Nakagami-m  and  Rice  PDFs  for  all  cases. 


Ionospheric  scintillation  using  GPS  over  South  America 


LISN  GPS  Network  &  3  more  Networks 


USN  GPS  Network  (in  green)  &  3  more  Networks  (in  red).  The  gray  arrows 
indicate  the  Jicamarca  and  Sao  Luis  digisondes,  the  diamond  indicate  the 
position  of  the  imager  at  Sao  Joao  do  Cariri. 


Ionospheric  scintillation  using  GPS  over  South  America 

Seasonal  percentage  of  scintillation  occurrence  for  3  stations 
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Estacao  Equatorial  de  Sao  Luis 
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Estacao  Equatorial  de  Jicamarca 


SET  OUT  NOV  DEZ 


Estacao  Equatorial  de  Christmas  Island 


JAN  FEV  MAR  ABR  MAI  JUN  JUL  AGO  SET  OUT  NOV  DEZ 


Black  histogram:  Observed  plumes  and  white  histogram:  Spread  F.  At  Sao  Luis  plume 
occurrence  is  from  September  to  March  with  a  peak  at  summer.  At  Jicamarca  the 
peaks  are  at  equinoces  and  for  Christmas  Island  is  on  August. 

The  difference  between  black  and  white  bars  represent  bottom-type  structures. 


Irregularity  prediction  model(Emanoel  Costa/Ricardo  Yvan) 


made  in  mapping  the  irregularities 
along  the  magnetic  field  lines. 


•  Radio  waves  passing  through  the  irregularities  diffract 
producing  signal  fading  and  strong  scintillations  even  at  L 
frequencies. 

•  The  contribution  in  the  model  was 


Irregularity  prediction  model(Emanoel  Costa/Ricardo  Yvan) 


Range-Time-Intensity  maps  from  the  Sao  Luis  coherent  scatter  radar  will  be  used  to 
predict  time  variations  of  the  scintillation  index  S4. 

The  model  use  the  s/n  ratio  measurements  by  the  VHF  radar  to  estimate  the  mean 
square  electron  density  fluctuation  within  the  corresponding  volume. 


We  are  using  zonal  drift  velocity  from  climatological  model  (Arruda  et  al.,  2006).  V  = 
150m/s  between  19-23  UT,  80  m/s  between  23-29  UT. 


Irregularity  prediction  model(Emanoel  Costa/Ricardo  Yvan) 

{AN1)  Mean  Square  Electron  Density 


S/N 


Signal  to  Noise  ratio 


Power  Spectral  Parameters 


Obtained  from  / 
C/NOFS  satellite  for 
2008 


Radar  Parameters 


Parameters 

Values 

L0[km] 

12.500 

0.080 

2.600 

q 

4.600 

si 

3.210 

Freq\M  H  z ] 

29.795 

Pt[kW] 

4.000 

Eff 

0.500 

dh\km\ 

2.500 

BW3dB[deg\ 

10.000 

Cain[dBi] 

26.000 

NFig\dB ] 

5.000 

Tsys[K} 

917.061 

Tsky[K\ 

10000.000 

BW[kHz] 

120.000 

2  70" 


Irregularity  prediction  model(Emanoel  Costa/Ricardo  Yvan) 
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Irregularity  prediction 

Ionospheric  irregularity  precursors  (prediction  tentative) 

Influence  of  MSTIDs  and  GWs  on  the  irregularity  generation  (Ricardo's  PhD) 
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SLZ  284/2009 


(a)  -  TEC  depletions 

(b)  -  TEC  disturbances(TIDs) 

(c)  -  01630  nm  images 

(d)  -  Bubble  separation 

(e)  -  Digisonde  altitude  and 
drift  velocity 

Wave  oscillations  of  about  30 
minutes  were  observeds  which 
give  evidency  of  gravity 
waves(GW) 


1 OO  200 


LOCAL  TIME 


Irregularity  prediction(Precursor) 


Irregular  base  layer  dynamics  before  the  bubble  triggering  using  radar  imaging 
(in  collaboration  with  Alam  Kherani  from  INPE) 
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With  time  this  layer  grows  upward  outside  tKe  shear 
region  and  becomes  a  bubble. 

The  bubble  can  grow  in  altitude,  detached 
from  the  base  layer. 


Irregularity  initiation  in  the 
base  layer  probably  generated 
by  the  Rayleigh-Taylor  and 
shear  instabilities 
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Irregularity  prediction(Precursor) 

Irregular  F  and  E  layers  dynamics  before  the  bubble  occurrence  using  VHF  Radar 

(a)  2005/09/29  (b)  2005/10/05 
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The  Doppler  velocity  inside  the  irregular  F  layer  shows  wavy  behavior. 

Similar  behavior  was  observed  in  the  irregular  E  layer  with  a  time  shift  prior  to  the  F  layer. 
And  this  E  layer  behavior  could  be  considered  as  a  precursor  of  the  bubble. 


Irregularity  prediction 


The  following  procedures  can  also  be  used  to  predict  irregularities: 

•  lonosonde  vertical  drift  calculations  to  establish  thresholds  to  trigger 
irregularities  (this  work  is  being  developed). 

•  Bubble  velocity  is  eastward  (during  quiet  period)  so  spaced  GPS 
receivers  /  VHF  receivers  in  the  zonal  direction  (SCINDA)  can  be 
employed  to  predict  bubble  occurrence  to  eastward  station. 


6  different  GNSS  receivers  campaign 

To  analyze  their  behavior  under  scintillation  condictions  (Septentrio, 

Novatel  4004B  and  GPS-Station6,  ASTRA, GEC-PLESSEY  Card  -  Cornell,  Stanford  -  U  Box  )  Keith  Groves,  Cesar 

Valladares,  Todd  Walter,  Geoff  Crowley,  Paul  Kintnert 


Campaing  of  6  different  GNSS  receivers(Stanford) 


6  different  GNSS  receivers  campaign 

Different  GNSS  receivers  behavior  under  scintillation  condictions  during  March  24  2013 

at  Sao  Jose  dos  Campos,  Sao  Paulo  -  Brazil 

Scintillation  monitors  -  Performance  comparison  (PRN1  at  March  24th) 


- U-Box 


22  22h30  23  23h30  24  24h30 

LT[Hours] 


Good  S4  agreement  for  the  6  receivers  and  for  moderate  scintillations  -  no  strong 
scintillations  tested 
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